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In this work we study a modified version of the two-dimensional Sznajd socio- 
physics model. In particular, we consider the effects of agents' reputations in the 
persuasion rules. In other words, a high-reputation group with a commom opinion 
may convince their neighbors with probability p, which induces an increase of the 
group's reputation. On the other hand, there is always a probability q = 1 — p of 
the neighbors to keep their opinions, which induces a decrease of the group's repu- 
tation. These rules describe a competition between groups with high reputation and 
hesitant agents, which makes the full-consensus states (with all spins pointing in one 
direction) more difficult to be reached. As consequences, the usual phase transition 
does not occur for p < p c ~ 0.69 and the system presents realistic democracy-like 
situations, where the majority of spins are aligned in a certain direction, for a wide 
range of parameters. 
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I. INTRODUCTION 

Social dynamics have been studied through statistical physics techniques in the last 
twenty years. Among the studied problems, we can cite models of cultural [1], language {2] 
and opinion dynamics 3-5] (for a recent review, see These models are interesting for 
physicists because they present order- disorder transitions, scaling and universality, among 
other typical features of physical systems [6]. 

Ising-type models, which can represent systems with binary variables and 1 (or ±1), 
have been used by physicists in many different areas, such as sociology, politics, marketing 
and finance J 
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. In 2000, an agent-based model was successfully applied to the dynamics of 
a social s^fcj . In this model, each agent carries one of two possible opinions, which may 
be represented by Ising spins. The aim of the Sznajd Model (SM) is the emergence of social 
collective (macroscopic) behavior due to the interactions among individuals, constituting 
the microscopic level of a social system jf|. 

The original SM consists of a chain of sites with periodic boundary conditions where 
each site (individual opinion) could have two possible states (opinions) represented in the 
model by Ising spins ("yes" or "no"). A pair of parallel spins on sites i and % + 1 forces 
its two neighbors, i — 1 and i + 2, to have the same orientation (opinion), while for an 
antiparallel pair + 1), the left neighbor (i — 1) takes the opinion of the spin i + 1 and 
the right neighbor (i + 2) takes the opinion of the spin i. In this first formulation of the 
SM two types of steady states are always reached: complete consensus (ferromagnetic state) 
or stalemate (anti- ferromagnetic state), in which every site has an opinion that is different 
from the opinion of its neighbors. However, the transient displays an interesting behavior, 
as pointed by Stauffer et al. [h]]. Defining the model on a square lattice, the authors in [uj] 
considered not a pair of neighbors, but a 2 x 2 plaquette with four neighbors. Considering 
that each plaquette with all spins parallel can convince all their eight neighbors, a phase 
transition was found for an initial density of up spins d = 1/2. 

A more realistic situation is to consider a probability of persuasion of each site. The 
Sznajd model is robust against this situation: if one convinces the neighbors only with some 
probability p, and leaves them unchanged with probability 1 — p, still a consensus is reached 



after a long time jg]. Models that consider many different opinions (using Potts' spins, 
for example) or defined on small-world networks were studied in order to represent better 



approximations of real communities' behavior (see |8j and references therein). In another 
work, in order to avoid full consensus in the system and makes the model more realistic, 
Schneider introduced opportunists and persons in opposition, that are not convinced by 
their neighbors [111, w hereas Stauffer considered that only individuals with similar opinions 
can be convinced [l2| . In a recent work, another formulation of the Sznajd model with 
limited persuasion was studied, considering that only high-reputation groups can convince 



neighbors 



In this work we analyze a modification of the Sznajd model with reputation studied in 



13J. We consider that the reputations may increase and decrease, depending if agents are 



persuaded or not by high-reputation groups. This fact makes the models' rules intrinsically 
dynamic, a common feature of evolutionary systems [14J |. The phase transitions as well as 
the relaxation times of the model are studied, and the relations between the quantities of 
interest and the parameters of the model are discussed. 



II. MODEL 



The model considered in this pap er is a modification of a recente work [l3j . Thus, 



we have considered the Sznajd model [10| defined on a square lattice with linear size L and 
periodic boundary conditions, where an integer number (R) labels each player and represents 
its reputation across the community. The reputation is introduced as a score which is time 
dependent. The agents start with a random distribution of the R values, and during the 
time evolution, the reputation of each agent changes according to its capacity of persuasion, 
following the model's rules. One time step in our model is defined by the following rules: 

1. We randomly choose a 2 x 2 plaquette of four neighbors. 

2. If not all four spins are parallel, leaves its eight neighbors unchanged. 

3. If all four center spins are parallel, we calculate the average reputation of the plaquette: 

! 4 
R = 4 Ri ' 

i=l 



where each term Ri represents the reputation of one of the plaquettes' agents 
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4. Then, we compare the average reputation R with the reputations of each one of the 8 
sites neighboring the plaquette. Thus, with probability p we follow the rule presented 



in [13|], i.e., given a plaquette neighbor j with reputation Rj, this individual j will 
follow the plaquette opinion if R > Rj. In this case, the reputation of each agent in 
the 2x2 plaquette is increased by 1 (in other words, the average reputation of the 
plaquette is increased by 1). 

5. On the other hand, with probability q = 1 — p the neighbor j keeps his opinion, even 
if we have R > Rj. In this case, the reputation of each agent in the 2x2 plaquette 
is decreased by 1. 

6. We repeat steps 4 and 5 for each one of the eight plaquette neighbors. 

Notice that agents belonging to a given plaquette may be rewarded or penalized in terms 
of reputation. A similar approach was analyzed in the Naming game, where a given agent's 
reputation is increased (decreased) depending on the outcome, i.e., success (failure) of his 
communication attempt with some other agent [3] . Agent-based models of animal behavior 



13, 



18J and human relations 19] have also considered reputation-like mechanisms in order 
to take into account different strengths or convincing powers. 

The above-mentioned rules define a competition between groups of high-reputation agents 
and hesitant individuals, that have a probability q = 1 — p of keeping their opinions despite 
the influence of plaquettes with high reputation. Due to above microscopic rules, one can 
expect difficulties to the system to reach the full-consensus steady state, with all spins 
oointing in a certain direction. Observe that for p = 1 we recover the model studied in 



131 ] . where democracy-like situations, with the majority of spins pointing in one direction, 
occur spontaneously, due to the dynamic rules considered in the model. Thus, concerning 
these realistic situations, we can expect similar results in our case, at least for values of the 
probability p near unity. 



III. NUMERICAL RESULTS 



In the simulations, the initial values of the agents' reputation follow a gaussian 
distribution centered at with standard deviation a = 5 [20|. We can start studying the 
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FIG. 1. Time evolution of the magnetization per spin for L = 73, some values of the initial 
density of up spins d and two different values of the probability p, namely p = 0.7 (left side) and 
p = 0.8 (right side). We can observe steady states with \m\ < 1, where there is no consensus in the 
population, as well as full-consensus states, where \m\ = 1. 
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where = ±1 and A = L is the total number of agents on the lattice. We exhibit in Fig. 
[T]the magnetization per spin as a function of the number of time steps for p = 0.7 (left side) 
and p = 0.8 (right side) and some values of the initial density of up spins d. Observe that 
the time needed for the system to reach the steady states increases for decreasing values of 
p. This result is easy to understand: when we decrease the value of the probability p, the 
probability q = 1 — p of the agents to keep their opinions increases, which makes stronger 
the competition between increasing and decreasing reputations. Observe that for p < 1 
situations with \m\ < 1 emerge spontaneously (see Fig. [[]), as in the p = 1 case [131 ] . On the 
other hand, full-consensus states with all spins pointing up or down are also possible in our 
model, as we can also see in Fig. [TJ Thus, the dynamics of competition between reputations 
leads to the emergence of democracy-like situations in our model. In addition, the results of 
Fig. [1] suggest that for p < 1 there are transitions between a phase where the full-consensus 
states with m = 1 are always reached and another phase where this consensus never occurs. 
We will analyze this phase transition in the following. 
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FIG. 2. Fraction / of samples which show all spins up when the initial density of up spins d is 
varied in the range 0.4 < d < 1.0, for L = 101 and typical values of the probability p. A phase 
transition can be already observed to occur at some critical point p c (0.68 < p c < 0.69), better 
analysed in the sequence. 

For the characterization of the phase transition, we have simulated the system for different 
lattice sizes L and we have measured the fraction of samples / which show all spins up when 
the initial density of up spins d is varied in the range 0.4 < d < 1.0. In other words, this 
quantity / gives us the probability that the population reaches consensus, for a given value 
of d. We have considered 1000 samples for L = 31 and 53, 500 samples for L = 73 and 101 
and 200 samples for L = 121. The numerical results for / as a function of d for L — 101 
and typical values of p are shown in Fig. |2j Observe that for values of p near 0.6 there is no 
transition, i.e., the consensus states are obtained only for d = 1 in these cases. For the formal 
characterization of the transition, we can study how the density values d corresponding to 
/ = 0.5 depends on the lattice size L. Taking the value p = 0.7, for example, we exhibit 
in Fig. [3] the results for / as a function of d for different lattice sizes L. We can estimate 
the values d(L, f = 0.5) from the crossing of the curves obtained from the simulations for 
each value of L and the dashed straight line in Fig. [31 Plotting these values as a function 
of L~ b one gets a straight line for b ~ 0.5 (see Fig. HJ left side). Fitting data, we obtained 
the equation 
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FIG. 3. Fraction / of samples which show all spins up when the initial density of up spins d is 
varied in the range 0.5 < d < 1.0, for p = 0.7 and some values of the lattice size L. The dashed 
line shows that the value of d corresponding to the level / = 0.5 increases for increasing values of 
L. 



d(L, f = 0.5) = 0.96007 - 1.3989 L" 5 , (2) 

which give us an estimate of the critical density d c ~ 0.96 in the thermodynamic limit 
L^ 1 — > 0. The scaling plot of the quantity / was obtained based on the equations 

f(d,L) = L- a f((d-d c )L b ) , (3) 
d c (L) =d c + cL- b . (4) 

The result is shown in Fig. H] (right side). The best collapse of data was obtained for 
b = 0.47 ± 0.03, a = 0.023 ± 0.002 and d c = 0.96007 ± 0.0031. 

Following the above- discussed process, we have estimated the critical points d c and the 
critical exponents a and b [see Eqs. ([2]), ([3]) and fll])] for other values of p. We have verified 
that the exponents a and b do not depend on p, considering the error bars. Thus, the 
criticallity of our model is denned by the critical exponents a ~ 0.02 and b ~ 0.5, and the 
phase transition occurs at specific critical points d c {p) which depends on the probability p. 
In addition, applying the process to identify the critical density d c [see Eq. (j2J)] for values 
smaller and near p = 0.7, we have found that d c {p = 0.69) = 0.9845 ± 0.0057, but for 
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FIG. 4. The value of d corresponding to the level / = 0.5 for different sizes L, d(L, f = 0.5), as 
a function of L~°- 5 (left side). Fitting data, we obtained d(L,f = 0.5) = 0.96007 - 1.3989 L -0 - 5 , 
which give us d c (p = 0.7) ~ 0.96 in the thermodynamic limit. It is also shown the scaling plot of 
/ (right side). The best collapse of data was obtained for a = 0.023 ± 0.002, b = 0.47 ± 0.03 and 
d c = 0.960 ± 0.003. 

p = 0.68 we have found that d c {p = 0.68) > 1, which suggests that for p < 0.69 the system 
doesn't undergo the usual phase transition of the Sznajd model. This defines a critical value 
p c ~ 0.69 below which the system does not present the phase transition between consensus 
and no-consensus states. Observe that this transition is robust with respect to the choice of 
different values of the standard deviation a of the initial distribution of agents' reputation 
(see Fig. [5]). Taking into account the critical densities d c (p) calculated for typical values of 
p, we exhibit in Fig. [6] the phase diagram of the model in the plane d c versus p. The points 
were estimated from the finite-size scaling analysis of the numerical results, whereas the line 
is just a guide to the eyes. Considering the thermodynamic limit L — > oo, the critical points 
d c separate a region where the full-consensus states with all spins up always occurs, for each 
realization of the dynamics (above the frontier), from a region where this kind of consensus 
never occurs (below the curve). 

We have also studied the relaxation times of the model, i.e., the time needed to find all 
the agents at the end having the same opinion. The distribution of the number of sweeps 
through the lattice, averaged over 10 4 samples, needed to reach the fixed point is shown in 
Fig. [7| (left side), for typical values of the probability p. We can see that the relaxation 
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FIG. 5. Fraction / of samples which show all spins up when the initial density d is varied in the 
range 0.4 < d < 1.0, for L = 23, p = 0.7, 1000 samples and typical values of the standard deviation 
a of the initial distribution of agents' reputation. This result shows that the increase of a does not 
change the behavior of /. 



time behavior is compatible with a log-normal distribution for all values of the probability 
p, which corresponds to a parabola in the log- log plot of Fig. [7J (left side). In addition, the 
average relaxation times increase for decreasing values of p, as a consequence of the above- 
mentioned competition of reputations. It is also exhibited in Fig. [7J (right side) the average 
relaxation time r [also over 10 4 samples, considering the relaxation times of Fig. [7J (left 
side)] versus latice size L in the log-log scale. We can verify a power-law relation between 
these quantities in the form r ~ L a , for large L, with a p— dependent exponent a. Observe 
that a increases for decreasing values of p, which is also a consequence of the competition 
of reputations. Notice that log-normal distributions for the relaxation times and power-law 



relations between r and L were also found in previous works on the SM 
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IV. FINAL REMARKS 

In this work we have studied a modified version of the Sznajd sociophysics model. In 
particular, we considered reputation, a mechanism that limits the capacity of persuasion of 
the agents. In other words, a high-reputation group with a commom opinion may convince 



10 



1| 1 ' r 

0.98- 

0.96- §-.._ 
0.94- 



T3 

0.92- 



0.9 no-consensus phase 

0.88- 
0.86 



consensus phase 



i; 



0.7 0.75 0.8 0.85 0.9 0.95 1 
P 



FIG. 6. Phase diagram of the model in the plane d c versus p. The points are numerical estimates 
of the critical densities d c , whereas the dashed line is just a guide to the eyes. In the region above 
the curve the system always reaches consensus with all spins up, whereas in the region below the 
curve this kind of consensus never occurs. 



their neighbors with probability p, which induces an increase of the group's reputation. On 
the other hand, there is always a probability q = 1 — p of the neighbors to keep their opinions, 
which induces a decrease of the group's reputation. These rules describe a competition 
between groups with high reputation and hesitant agents, which makes the full-consensus 
states (with all spins pointing in one direction) more difficult to be reached. The consequence 
is the occurrence of democracy-like situations, where the majority of spins point in one 
direction, for a wide range or parameters. 

We have observed a log-normal-like distribution of the relaxation time, i.e., the time 
needed to find all the agents eventually having the same opinion. In addition, the average 
relaxation times grow with the linear dimension of the lattice in the form r ~ L a , where a 
depends on the probability p. We have identified the critical density d c (p) for some values of 
p by the finite-size scaling analysis of the fraction / of samples which show all spins up when 
the initial density of up spins d is varied. Our numerical results suggest that the system 
does not undergo the usual phase transition for p < p c ~ 0.69. 

It is important to notice that the inclusion of the new ingredient, i.e., reputation, makes 
the Sznajd model intrinsically dynamic. Whereas in the previous formulations of the SM 
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FIG. 7. Log-log plot of the histogram of relaxation times for L = 53, d = 0.9 and some values of 
p, obtained from 10 4 samples (left side). The observed parabolas in the log- log plot indicate log- 
normal-like distributions for all values of the probability p. It is also shown the average relaxation 
times r as functions of the lattice size L (right side). Fitting data, we obtained the power-law 
behavior r ~ L a , where the numerical values of the exponents a are exhibited in the figure. 

the microscopic rules are static in a way that the dynamics is only related with the evolution 
of the system, in our case the SM becomes completely dynamic since the rules evolve with 
the system, i.e., the rules change according to the succession of visited states. 
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